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We study the formation of dark-bright vector soliton pairs in nonlocal Kerr-type nonlinear medium. We 
show, by analytical analysis and direct numerical calculation, that in addition to stabilize vector soliton pairs 
nonlocal nonlinearity also helps to reduce the threshold power for forming a guided bright soliton. With help of 
the nonlocality, it is expected that the observation of dark-bright vector soliton pairs in experiments becomes 
more workable. © 2008 Optical Society of America 
OGIS codes: 190.5530, 190.3270. 



Recently the study of nonlocal nonlinearity brings new 
features in solitons.^ such as modification of modulation 
instability^ and azimuthal instability,^ suppression of 
collapse in multidimensional solitons,'* change of the soli- 
ton interaction,^ and formation of soliton bound states.^ 
Nonlocal effect comes to play an important role as the 
characteristic response function of the medium is compa- 
rable to the transverse content of the wave packet. Ex- 
perimental observations of nonlocal response also have 
been demonstrated in various systems, such as pho- 
torefractive crystals,^ nematic liquid crystals,^ thermo- 
optical materials,^ and ^^Cr Bose-Einstein condensates 
with strong dipole-dipole interaction.^*' 

For nonlinear local media, dark-dark, bright-bright, or 
dark-bright soliton pairs can exist in normal or anoma- 
lous dispersion region in vector settings in contrast to 
the scalar models. ^^'^^ Especially with the help of a 
dark soliton, a bright soliton in normal dispersive me- 
dia is formed through soliton-induced guiding effect. 
With nonlocal nonlinearity, a large number of multi- 
hump multi-component vector solitons are found with a 
remarkable stabilization.^'* For nonlocal Kerr-type non- 
linear medium, the nonlocality is known to improve the 
stabilization of solitons due to the diffusion mechanism 
of the nonlinearity. The price to pay is that nonlocal 
solitons also need to increase their formation power to 
compensate the diffusion effect in nonlocal materials. In 
this Letter, we study vector solitons with dark-bright 
pairs in nonlocal nonlinear media. We reveal that stable 
bright solitons can be formed in normal dispersive re- 
gion guided by dark soliton backgrounds. Moreover we 
find that the nonlocal nonlinearity also helps to reduce 
the threshold power for such a guided bright soliton due 
to the combination of nonlocality and vectorial coupling. 

We consider two mutually incoherent wave packet 
propagating along the ^ axis within a nonlocal Kerr- 
type nonlinear medium. The governing equations of the 
vectorial Manakov system which consists of two vector 



components U and V are given by, 

.dU Id^U ,^ ^ 
i—--—+nii,7^)U = 0, 
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+ n(^,77)F = 0, 



(1) 
(2) 



2 drj^ 

n{^,v)= r Riv-v')m' + \V\')drf, (3) 

(4) 



where rj is the transverse coordinates, n{^,r]) is the re- 
fractive index profile induced by the exponential-type 
diffusion kernel function R{r]) responding to the intensity 
soliton intensity.*^ The coefficient d stands for the de- 
gree of nonlocality which governs the diffusion strength 
of refractive index. With the nonlocal vector model in 
Eq. (HJIl]), stationary solutions in form of V{£_,ri) = 
u(77)exp(i^„^) and U{£,,ri) = u{ri)exp{i^u£,) are assumed 
to be the solution of dark-bright vector soliton pairs with 
real propagation constants, /i„ and respectively. The 
two component solutions of dark-bright vector soliton 
pairs are subject to the boundary condition u{±oo) = 
and f (±cx)) = ±y/JI^. 

The solutions of dark-bright vector soliton pairs in 
local and nonlocal media, as well as the refractive in- 
dex profiles, are shown in Figure [T] (a) and (b). The 
dependence of the power for bright component, defined 
in Eq.Q, and its propagation constant is plotted in 
Fig. [1] (c). As the case in scalar model, a bright soli- 
ton in nonlocal media has higher cutoff potential than 
the local one due to the diffusion of the nonlinear in- 
dex. In Figure [21 the relations between the propagation 
constant for bright soliton /i„ and the degree of nonlo- 
cality d at different fixed powers are shown. It can be 
seen that the propagation constant of the bright one in 
a dark-bright soliton system is growing as the degree of 
nonlocality increases, for the bright component sees a 
potential directly from the nonlinear index that is raised 
due to the diffusion nonlocality in response to the inten- 
sity sum of the dark-bright soliton pair. But in contrast 
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to the scalar soliton in local mcdia/^ the bright soliton 
guided by a dark soliton in vector model requires lower 
forming power in nonlocal region as the same propaga- 
tion constant is concerned, as the marked points A and 
B shown in Fig. [1] (c). 

In addition to direct numerical treatments, we exam- 
ine the formation power of bright soliton in nonlocal 
media analytically by variational methods. To simplify 
the analysis, a constant dark pulse is assumed as the 
background and the Lagrangian equation for the bright 
soliton in a vectorial nonlinear nonlocal system has the 
form, 



L 



uiu) + l\u,f + l-\uf + \uf\v\' 



M + l\u\'{u,,u* + u;^ u) + If/P^lv'P 



(5) 
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Fig. 2. Relations between the propagation constant for 
bright soliton and the degree of nonlocality d at dif- 
ferent fixed powers, P = 0.5 (solid line) and P = 1.0 
(dashed line). 



where the subscriptions ^ and 77 stand for derivative with 
respect to longitudinal and transverse coordinates. The 
terms within the brackets multiplied by d in the sec- 
ond line in Eq. ([5]) represent the nonlocal index response 
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Fig. 1. Intensity profiles of bright soliton (solid line), 
dark soliton (dashed line) , and the refractive index (bold- 
dashed line) for the local (a) and nonlocal media (b) at 
the points A and B in (c), respectively. The bifurcation 
curves of fundamental dark-bright soliton pairs in local 
{d = 0, dashed line) and nonlocal media (d = 1, solid 
line) are shown in (c), where points A and B arc marked 
with /x„ = 0.8 and /z^ = 1. 



which can be calculated through following expansion 
1 , <9™|i7|2 



— ' m! drf 

ni—O ' 

\U\^ + \V\^+d 



(6) 



dr]' 



where 



dcj™ 



is the expansion of the Fourier transform, H{lo), of the 
kernel function R{ri) in Eq. To solve the Lagrangian 
equation in Eq. ([S]), we use following solution ansatz for 
the bright and dark solitons, 

U{r]) = Ausech{ri/au)exp{i4)u + icufl^), (7) 
V{r]) = Aytanh{ri/av)exp{i(j)v + iCyT]'^), (8) 



where the parameters Aj , aj , and Cj , (j = U, V) 
are amplitude, width, phase, and chirp for bright and 
dark solitons, separately. By assuming that dark soliton 
is invariant to the change of the degree of nonlocality in 
the low nonlinear limit, a set of Euler-Lagrangian equa- 
tions for Au, au,4>u and Cj can be obtained. Furthermore, 
we assume that in steady state, (j)u = /x^ is a constant, 
and Cj is zero for chirpless soliton solutions. Then for a 
set of propagation constants of dark-bright soliton pairs, 
/i^ and fj,u, an approximate linear dependence of bright 
soliton power versus nonlocality is derived as. 



P= |C/|^d?7 = 
V2/it, - 2/i„ 



2A^au, 



0.13 



0.877 



(9) 



1.57 



In this first-order approximation, Eq. Q indicates that 
the forming power of bright soliton guided by a dark 



2 



Numerical 

Theoretical 




Fig. 3. Threshold power of bright sohtons versus the de- 
gree of nonlocahty for ^„ — 0.8 and fiy = 1.0. SoUd and 
dashed Hnes are calculated by numerical and variational 
methods, respectively. 

soliton in vector model decreases as the degree of the 
nonlocahty increases. In addition, when d = 0, Eq. ^ 
reduces to the case of soliton solutions in the local me- 
dia. In Fig. [HI we show the dependence of threshold 
power for forming bright solitons with the degree of non- 
locality both by direct numerical simulation of Eq. ([l]|4]) 
and the Lagrangian equation in Eq. ([5]), which is con- 
sistent with the results in Fig. [ijc). In this case we fix 
fly to 1 since it is associated with the dark component 
which is given by the boundary conditions and can be 
scaled out. Conceptually, as the degree of nonlocahty in- 
creases, the tendency for refractive index to advance to 
the region of lower light intensity grows stronger. Even 
though the dark pulse almost remains unchanged, the 
existence of bright pulse drives out index flow. Conse- 
quently the index modulation induced by the soliton pair 
becomes shallower, and the nonlinearity required to form 
the bright soliton decreases. The power of bright soliton 
decreases in a dynamical balance with the refractive in- 
dex flow. This implies that the threshold power to form 
a bright soliton guided in the dark background can be 
reduced with nonlocal interaction. 

The stability of dark-bright soliton pairs in nonlocal 
nonlinear media is analyzed by standard linear stability 
analysis with introduction of perturbations solutions, i.e. 
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where the small perturbations grow at the rate of the 
real part of A. As the case in local media, dark-bright 
soliton pairs in Manakov model are stable in the nonlocal 
nonlinear media. 

In conclusion, we study the formation of dark-bright 
soliton pairs in vectorial nonlocal nonlinear model ana- 
lytically and numerically. We find that in addition to the 
stabilization of vector soliton pairs, nonlocal nonlinear- 
ity also helps to reduce the threshold power for forming 
a guided bright soliton due to the dynamical balance be- 
tween the nonlinearity and nonlocal induced refractive 



index flow. With a constant background of dark soli- 
tons, our analytical model shows a linear dependence of 
the formation power for bright solitons on the degree 
of nonlocahty and also matches the numerical simula- 
tions very well. With the reduction of forming threshold 
power, we believe that our results arc very useful for the 
observation of dark-bright vector soliton pairs in nonlo- 
cal nonlinear media. 

Authors are indebted to Yu. S. Kivshar, W. 
Krolikowski, O. Bang, and A. S. Desyatnikov for useful 
discussions. This work is supported by the National Sci- 
ence Council of Taiwan with the contrast number NSC- 
95-2120-M-001-006. 
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